DOC  FILE  COPY  AD  A0938  7  3 


AFOSR-TR-  8  0*134 


FINAL  SCIENTIFIC  R.1P0RT 


SUBJECTED  TO  INTENSE  ELECTRON  BOMBARDMENT 


GRANTS  AFOSR-80-0035 


Prepared  for  the 
Air  Force  Office  of  Scientific  Research 
by 

'"^Arizona  Research  Labs,;  I  ^ 
Professor  Marian  0.  Scully,  Principal  Investigator 


November,  19S0 


///ot  i  33 


Approved  for  publio  release 
distribution  unlimited. 


/{REPORT  DOCUMENTATION  paqe 


RUAD  INST KUCTI ON'S 
HKI-OKL  COMI'M-.T 'I NT.  1-ORM 


10  1  o.’  Scully,  J.  K.^McIvcr/ G.  T.yjloore 


/AFOSR- 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS  ICE  7-7.* GRAM  LLr.HCNT,  PROJECT,  TASK 

AREA  it  A’ORr;  UNIT  NUMBERS 

University  of  Arizona  _ 

Tucson,  Arizona  85721  YjjJ  S’'jz?.2  ///^: 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS  #  ^  •  REPORT  DAT  _ 

//V^r  (if  |  No 

,  ~  _  ...  „  V>»—  ^I3l  NUlJBCP'Or  PACES 

/)iu ■>!/;-  ^  ,  0)0  29 

Ut~'MONITORING  AGEN.CY  NAME  »  AODRESST//  dillcrent  Irom  Controlling  Ollirc)  IS.  SECURITY  CLASS,  (ol  this  reporlj 


0//OJO- 

‘"V;  REPORT  DAT  £.  .  ,w  . 

//VyF  [  If  |  No 

V  p.131  nUnIBER  OE  P  AGES 

29 


Jl- 


m  a?- .  (' 


S/s/tiO 


1  5«.  DECL  ASSl  F|C  ATION'  DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (of  thia  Report) 


Approved  for  public  release  ; 

distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  in  Block  20,  It  different  from  Report) 


18.  SUPPLEMENTARY  NOTES 


M9.  KEY  fO°D£  { Continue  on  reverse  .tide  if  necessary  and  identify  by  block  number) 


20  ABSTRACT  f Continue  on  reverse  side  if  necessary  and  identify  by  block  mimber) 

Several  aspects  of  the  generation  of  x-rays  in  a  crystal  were  investigated  in 
order  to  explain  some  of  the  experimental  results  observed  by  Das  Gupta.  Sequen¬ 
tial  excitation  of  planes  of  two-level  atoms  in  a  crystal  was  shown  .to  produce 
stimulated  line  radiation  whose  intensity  varied  exponentially  with  the  incident 
electron  beam  current.  This  result  agreed  qualitatively  with  the  experimental 
da' a.  Swept  excitation  was  also  used  to  explain  the  increase  in  the  amount  of 
line  radiation  produced  at  a  given  election  beam  current  when  the  election's 
(Continued  on  back/  t  I 

DD  1473  J3  -Ufa-  .......  A 


(  //-  /s  y) 


mm 
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energy  was  changed  from  8.95  KeV  to  16.997  KeV.  The  difference  between  the 
lifetime  of  an  excited  two-level  atom  in  free  space  and  the  lifetime  of  the 
same  atom  in  a  lattice  of  nonresonant  two-level  atoms  was  also  investigated 
The  resulting  difference  was  attributed  to  Bragg  scattering  of  the  emitted 
radiation  in  the  crystal  and  the  bulk  properties  of  the  material. 
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ABSTRACT 

Several  aspects  of  the  generation  of  x-rays  in  a  crystal  were 
investigated  in  order  to  explain  some  of  the  experimental  results  observed 
by  Das  Gupta.  Sequential  excitation  of  planes  of  two-level  atoms  in  a 
crystal  was  shown  to  produce  stimulated  line  radiation  whose  intensity 
varied  exponentially  with  the  incident  electron  beam  current.  This 
result  agreed  qualitatively  with  the  experimental  data.  Swept  excit.  n 
was  also  used  to  explain  the  increase  in  the  amount  of  line  radiation 
produced  at  a  given  electron  beam  current  when  the  electron's  energy 
was  changed  from  8. 95  KeV  to  16.997  KeV.  The  difference  between  the 
lifetime  of  an  excited  two- level  atom  in  free  space  and  the  lifetime  of 
the  same  atom  in  a  lattice  of  nonresonant  two-level  atoms  was  also 
investigated.  The  resulting  difference  was  attributed  to  Bragg  scattering 
of  the  emitted  radiation  in  the  crystal  and  the  bulk  properties  of  the 
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The  purpose  of  this  research  has  been  to  explain  the  experimental 
results  obtained  by  Das  Gupta.1  These  results  include  a  nonlinear  increase 
in  the  amount  of  radiation  emitted  from  a  copper  target  and  a  decrease 
in  the  Kq  linewidth  when  the  intensity  of  the  incident  electron  or  x-ray 
beam  is  increased. 

Our  research  has  been  oriented  towards  two  major  areas  of  investigation. 
In  the  first,  sequential  excitation  of  the  atoms  in  a  crystal  has  been 
invoked  as  a  possible  means  of  explaining  the  nonlinear  increase  in  the  Kq 
intensity.  The  results  of  calculations  for  a  crystal  that  consists  of  a 
linear  array  of  atoms  shows  that  the  intensity  of  the  emitted  radiation 
increases  exponentially  as  the  intensity  of  the  pump  beam  is  increased.  This 
result  is  in  qualitative  agreement  with  the  experimental  data.  Quantitative 
comparisons  are  meaningless  until  further  work  is  completed.  In  particular, 
the  effects  of  the  three-dimensional  crystal  lattice  on  the  direction  and 
linewidth  of  the  emitted  radiation  must  be  taken  into  account.  In  order 
to  study  the  linewidth,  the  second  area  of  research  was  directed  towards  an 
investigation  of  the  properties  of  the  radiation  spontaneously  emitted 
by  an  excited  atom  in  a  crystal.  The  preliminary  results  of  this  calcula¬ 
tion  show  that  the  lifetime  of  an  excited  two-level  atom  in  a  crystal  of 
nonresonant  two-level  atoms  is  different  from  the  lifetime  of  the  same 
excited  atom  in  free  space.  This  work  may  have  also  provided  the  first 
steps  towards  understanding  the  linewidth  narrowing  observed  by  Das  Gupta.1 
Further  work  on  the  stimulated  emission  linewidth  of  an  excited  atom  in 
a  crystal  is  necessary  before  the  relevance  of  this  approach  can  be 


determined. 
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The  nonlinear  rise  in  the  amount  of  K  and  K  radiation  emitted 

cq  a2 

1  2 

from  a  copper  target  has  been  seen  using  both  electrons  and  x-rays 
to  excite  the  copper  atom.  Because  beams  of  both  charged  and  neutral 
particles  produce  the  same  nonlinear  response,  the  mechanism  responsible  for 
this  increase  is  to  first  order  independent  of  the  method  of  excitation. 

In  particular,  this  means  that  those  effects  that  rely  exclusively  on  a 
charged  excitation  beam  need  not  be  considered  at  first.  The  elimination 
of  charged  particle  effects  such  as  stimulated  quasi-Cherenkov  radiation 
and  channeling  greatly  simplifies  the  problem. 

radiation  is  produced  in  a  two  step  process.  In  the  first  step  a  Is 
electron  is  removed  from  the  atom  by  the  pump.  A  2p3^2  electron  then  makes  a 
radiative  transition  to  a  Is  state  producing  a  photon.  The  dominant 
competing  process  is  Auger  recombination  where  the  energy  released  by  the 
2p3/2  -*■  Is  transition  removes  a  higher  level  electron  from  the  atom  thus 
leaving  the  atom  doubly  ionized.  The  Auger  recombination  time  t  is  'vlO  15 

a 

sec.  This  means  that  for  an  initially  excited  medium  only  a  volume  of 

dimensions  cra  'v  3x10  5cm  could  participate  in  stimulated  scattering.  This 

is  much  smaller  than  the  characteristic  dimension  of  a  copper  crystal 
_ ,  2 

L'V'IO  3cm  reported  by  Das  Gupta.  The  limit  imposed  by  the  Auger  recombination 
can  be  removed  if  the  pumping  is  considered  as  swept  gain  excitation.  This 
is  the  model  of  the  excitation  mechanism  that  we  have  chosen  for  the  Das 
Gupta  experiment.  The  experiment  can  then  be  modeled  in  terms  of  a  two 
level  system. 
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Figure  1 

where  2  is  the  state  with  an  electron  missing  in  the  Is  shell,  1  is  the 
state  with  a  filled  Is  shell;  y  is  the  Auger  recombination  rate;  Y2>  Yi 

a 

are  the  decay  rates  of  the  two  states;  and  A  is  the  rate  at  which  state  2 
is  populated. 

The  density  matrix  equation  of  motion  along  with  Maxwell's 
equations  form  a  closed  set  of  equations  that  describe  this  two  level 
system.  These  equations  can  be  solved  in  the  weak  field  limit  (see  Appendix 
A).  The  resulting  expression  for  the  intensity  of  the  radiation  predicts 
an  exponential  increase  as  a  function  of  A.  This  is  at  least  in  qualitative 
agreement  with  Das  Gupta's  results.  (If  Das  Gupta's  data  on  the  variation 
of  the  intensity  as  a  function  of  pump  intensity  is  plotted  on  semilog 
graph  paper  then  the  resulting  points  very  nearly  lie  on  a  straight  line. 

This  is  to  be  contrasted  with  Das  Gupta's  plot  ).  A  quantitative  comparision 
of  the  predictions  of  the  model  and  the  experimental  data  has  not  yet  been 


made. 


Das  Gupta  has  recently  told  us  that  if  the  copper  sample  is 


excited  with  16.997  KeV  x-rays  instead  of  8.95  KeV  x-rays  then  the  nonlinear 

response  of  the  Cu  K  line  intensity  to  the  pump  beam  intensity  is  even  more 

pronounced.  The  two  energies  8.95  KeV  and  16.997  KeV  for  the  x-ray  pump 

beam  corresponds  to  E.  =  8.95  KeV  which  is  the  ionisation  energy  for  a  Is 
1  ion 

electron  and  E  =  16.997  KeV  =  E.  +  E  where  E„  =  8.047  KeV  is 

ion  i\  iv 

C»1  dj 
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the  energy  of  a  Cu  K  photon.  Because  the  pump  beam  is  incoherent,  the 
experiment  with  E  =  16.997  KeV  can  be  described  in  terms  of  two  decoupled 
two- level  systems  as  follows: 


E(z,t,w2-o>i) 


Figure  2 

where  =  16.997  KeV,  iiCwg  -  o>i)  =11(0)4  -  W3)  =  8.047  KeV  and-h(o)3  -  wj) 

8.950  KeV.  The  two  level  system  comprising  states  1  and  2  is  the  same  as 
described  in  Fig.  1.  The  remaining  two  level  system  is  composed  of  states 
4  and  3  .  State  4  is  an  energy  state  in  the  continuum  16.997  KeV  above  the 
Is  state  of  a  copper  atom.  It  is  pumped  at  the  same  rate  as  state  2  because 
in  the  experiment  the  removal  of  a  Is  electron  by  a  16.997  KeV  photon  not  only 
creates  an  atom  in  state  2  but  also  leaves  the  ejected  electron  with  a  kinetic 
energy  of  8.047  KeV  in  state  4  .  State  3  is  an  energy  state  at  the 
ionization  threshold.  The  density  matrix  equation  for  the  two  two-level 
systems  plus  Maxwell's  equations  again  form  a  complete  set.  The  resulting 
expression  for  the  intensity  of  the  copper  K  radiation  again  predicts  an 
exponential  increase  as  a  function  of  the  pumping  rate  (see  Appendix  B)  in 
the  weak  field  limit.  For  the  case  where  there  is  no  broadening  mechanism 
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present ,  the  difference  between  the  gain  coefficient  when  the  energy  of  the 
pump  beam  is  16.907  KeV  and  the  gain  coefficient  when  the  energy  of  the 
pump  beam  is  8.950  KeV  is  given  by 

O 

u>2 1  Pi 2^ 

a(16.997  KeV)  -  a(8.950  KeV)  =  — —  6 (?) f (C)  , 

where  pj 2  is  the  dipole  moment,  N  is  the  number  of  atoms  per  unit  volume, 

5  =  t  -  z/c,  0  is  the  unit  step  function  and  f(£)  =  (e"Y|+3  ^  -  e‘Yl+  ?)/(y2  -  yj). 
The  predictions  of  the  above  expression  have  not  yet  been  compared  with  the 
expei iment  because  Das  Gupta  has  not  yet  completed  the  experiment. 

In  order  to  make  a  quantitative  comparison  between  our  two  models  and 
the  results  of  the  experiment  it  is  necessary  to  understand  what  effect  the 
crystal  has  on  an  excited  atom.  P.  S.  bee  and  Y.  C.  Lee^  have  already  shown 
that  if  an  excited  two-level  atom  is  placed  in  a  crystal  consisting  of  identical 
resonant,  two- level  atoms  then  the  spontaneous  decay  rate  of  the  excited 
state  is  vastly  different  from  its  free-space  value.  For  the  case  of  the 
Das  Gupta  experiments,  the  K^  radiation  emitted  by  an  initially  excited  atom 
is  not  resonant  with  the  unexcited  copper  atoms  in  this  crystal.  A  model 
relevant  to  the  conditions  of  this  experiment  consists  of  a  single  excited 
two-level  atom  in  a  crystalline  array  of  nonresonant  two-level  atoms. 

Approximate  solutions  for  the  probability  amplitudes  of  the  various  states 
of  the  system  have  been  found  with  the  aid  of  the  multiple  time  scale 
perturbation  theory  of  Lindstedt  and  Poincare.^  With  this  method  the 
probability  amplitudes  have  been  solved  to  fourth  order  in  the  coupling 
constants  (see  Appendix  C  for  details).  Of  particular  interest  is  the 
complex  frequency  of  the  amplitude  for  the  initially  excited  atom.  The 
expression  for  this  frequency  consists  of  five  terms.  Three  of  these  terms 
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correspond  to  the  complex  frequency  that  would  exist  if  the  atom  were  in 
free  space.  The  remaining  two  terms  take  into  account  the  fact  that  the 
atom  is  in  a  medium.  One  of  these  terms  is  independent  of  the  structure 
of  the  material  whereas  the  other  takes  into  account,  the  lattice  structure. 
It  is  this  latter  term  which  is  of  interest.  It  takes  into  account  the 
Bragg  scattering  of  the  emitted  radiation  by  the  crystal.  A  rough  estimate 
of  the  imaginary  part  of  the  complex  frequency  indicates  it  is  less  than 
1%  of  the  Weisskopf-Wigncr  linewidth  and  is  therefore  not  important  in 
most  cases.  The  effect  of  the  crystal  on  the  stimulated  linewidth  has 
yet  to  be  determined. 

A  report  of  our  work  on  spontaneous  emission  of  an  excited  atom  in  a 
crystalline  array  of  nonresonant  identical  atoms  was  recently  presented 
at  the  International  Quantum  Electronics  Conference  in  June,  1980. 

Our  present  results  explain  qualitatively  some  of  the  results  of 

Das  Gupta's  experiments.  Quantitative  comparisons  between  the  prediction 

of  this  model  and  the  experiments  cannot  be  made  until  the  effects  of  the 

crystal  are  built  into  the  swept  gain  model.  In  particular,  the  Bragg 

scattering  of  emitted  x-ray  radiation  must  be  taken  into  account  in  order 

to  make  an  accurate  prediction  of  the  K  radiation  measured  in  the 

a 

experiment.  This  objection  does  not  hold  for  the  prediction  of  the  increase 
in  the  gain  coefficient  resulting  from  increasing  the  energy  of  the  pump 
beam  from  8.95  KeV  to  16.995  KeV.  A  quantitative  comparison  of  our 
prediction  and  the  results  of  the  experiment  await  the  completion  of  the 
experiment. 
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APPENDIX  A 

In  the  experiment  an  electron  or  an  x-ray  photon  ionizes  a  copper  atom 
by  removing  an  electron  from  the  Is  shell. *  The  vacancy  in  the  Is  shell  is 
removed  either  by  a  radiative  transition  of  an  electron  from  a  higher  level 
(2p  electron  in  the  case  of  Ka  radiation)  or  by  Auger  recombinations.  The 
Auger  recombination  is  not  a  radiative  transition.  The  energy  released 
by  a  2p  electron  going  to  the  Is  vacancy  is  used  to  remove  a  higher  lying 
electron  from  the  atom  thus  leaving  it  doubly  ionized.  Since  the  electrons 
involved  in  these  two  processes  are  either  in  the  Is  or  2p  shell  of  the 
copper  atoms,  the  atom  can  to  a  good  approximation  be  modeled  as  a  two- 
level  atom. 

The  characteristic  life  time  of  the  excited  atom  before  it  decays  by 
Auger  recombination  is  f  =  10  15  sec.  This  means  that  stimulated  emission 

a 

for  a  group  of  excited  atoms  will  only  occur  for  those  atoms  that  are 

within  ct  'v  10  5cm  of  each  other.  This  is  much  smaller  than  the  characteristic 

a 

size  10  3cm  of  a  crystal  in  the  sample  used  in  this  experiment. 

One  method  of  increasing  the  distance  over  which  stimulated  emission 
can  occur  is  to  sequentially  excite  the  atoms  along  some  axis.  This  is 
what  is  referred  to  as  swept  gain  amplification. 

The  various  processes  that  occur  in  the  tiro- level  atom  can  then  be 
represented  by  the  diagram: 
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where  state  2  is  the  excited  state  of  the  atom  (an  electron  is  removed  from 
the  inner  shell),  state  1  is  the  final  state  of  the  atom  (Is  shell  is  filled), 
Y  is  the  Auger  decay  rate,  y  and  are  the  decay  rate  for  the  lower  and 
upper  states,  E(z,t)  is  the  incident  electromagnetic  field  of  frequency  v 
and  A  is  the  rate  at  which  the  atoms  are  excited.  Because  these  atoms  are 
in  a  crystal  A  can  be  a  function  of  z. 

This  system  is  described  by  the  equations  of  motion  for  the  density 
matrix  elements 

p22  =  A  (z)S  (t  -  |-)  -  (y  +Y2)  p  22-  r  CV2  lP  12  ■  v12P2l),  (A.l) 


Pll  =  Y  P22  ~  YlPll  "  }~  (V12P2I  "  V2lPl2), 


p2l  =  -  i(w2  +  f)p21  +  r  (P22  "  Pll)V21 


(A.  2) 


(A.  3) 


and  the  wave  equation 


32  ...  1  a2  v,  _  a2P 

E(z,t)  -  E(z,t)  Pq  > 

az2  c2  at2  at2 


P(z,t)  =  N)j2  1  (p  2 1  +  P12). 


(A. 4) 


(A.  5) 


where  V21  =  -  y21E(z,t),  p21  =  e<2 1 x 1 1> ,  T  =(yj+  Ya  +  Y2)/2,  N  is  the  number  of  atoms 
and  fiu>2i  is  the  difference  in  energy  between  states  2  and  1. 

The  electromagnetic  field  and  polarization  field  can  be  rewritten  in 
terms  of  a  rapidly  oscillating  part  and  a  slowly  varying  part  as 


E(z,t)  =  £(z,t)  cos(kz  -  vt)  , 


(A.  6) 
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P(z,t)  =  S(z,t)  sin(kz  -  vt)  ,  (A. 7) 

with  the  conditions 

,  _  3£  aa  .  ^  9S  3S 

k&  »  ^  V  »  ks  »  ^  VS  »  ^  .  (A.8) 

We  shall  assume  exact  resonance  =  v  from  now  on.  The  substitution 
of  Eqs.  (A. 6)  -  (A.8)  into  Eq.  (A. 4)  yields 


where  £  =  t  -  —  and  G(C)  =  1  when  £>0  and  =  0  when  £>0.  The  combination 
of  Eqs.  (A. 9)  -  (A. 11)  and  (A. 12)  then  give  an  equation  for  the 
amplitude  of  the  electromagnetic  field 


H  +  c  H  =  2r  fdt'0(C')£(z,t')e-r(t't') 


f  -yi£'  -(y  +Y2K1 

-Cy  +Y2)?'  Ya[e  -e 

X  e  +  - - r - 

Yl  -  (Ya  +  Y2) 


(A. 13) 


where  r =  vpfiNA (z)/ce0h.  This  equation  can  be  rewritten  in  terms  of  the 


amplitudes  squared 


(§7  +  Flt)£2(z>t)  =  r(z)£Cz*t)  )tdt'  e  (Ya+Y2)E 

[  “YlC '  -(Ya+Y2K'i  '  0 

Ve  -e  -rct-t') 

+  -  tT-  CYa*Y2)  JK'VW  (  D 


!-(y  +Y2)5f  Y_ 
e  1  +  — ~ 

Y2+Y, 


/  _ 


t  (  Ya*Y2-Yl)T 


Ye  rt, 

*  £(z  t-x)  +  - - r-  die 

Yl“(Ya+Y2) 


Yl"(Ya+Y2) 

* - 2 - T 


£(z,t-x)l  . 


(A. 14) 


Because  and  Ya*  are  small  in  comparison  to  the  time  for  d?(z,t)  to 
change  significantly  (the  medium  is  excited  by  an  essential  continuous  beam) 
then  <?(z,t-x)  =  <?(z,t)  so  that  Eq.  (14)  becomes 


(infifx)82'2-*’  ■ 


F(0  =0(4)  1  + 


Ya+Y2-Yl 


CYa+Y2-Yl)T/2 
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Y1-(Y+Y2)  J 


Yl  -  Cy„+Y2  ) 

v  a 

f£  - ^ -  t 

dre 


If  we  change  variables  from  (z,t)  to  (5  =  (t  -  — ) ,  z '  =  z)  then 


3_  =  3__  3Z  +  3 3zJ_  _  3 _1  3_ 

3z  3^  3z  3z '  32  3z'  '  c  3?  * 


3_  _  3_  3£  _  3_ 
3t  ~  35  3t  35 


so  that 


L  +  I L  + 

32  c  3t  3z ' 


and 


Ip-  £2(z',5)  =  r(z '  )F(5)S2  (z '  ,£)  , 


which  has  a  solution 


£2rz‘,5) 


S2(0,5)  exp{ 


z ' 

r(x)dxF(5)}  . 


0 


(A. 


(A. 
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If  A(z)  =  Af(z)  then  we  see  that  Eq.  (A. 17)  predicts  that  the  intensity 

of  the  radiation  will  increase  exponentially  as  the  pumping  rate  A  is 

increased.  This  result  is  in  qualitative  agreement  with  the  results  of 

1 

Das  Gupta's  experiment. 

It  is  important  to  note  that  Das  Gupta  has  fitted  his  data  to  a  poly¬ 
nomial  in  the  incident  beam  intensity.1  An  equally  good  fit  can  be 
obtained  by  fitting  the  data  to  an  exponential  in  the  incident  beam  intensity. 
This  is  particularly  true  for  the  results  of  the  data. 

Although  this  model  qualitatively  agrees  with  the  experiment,  a  quantita¬ 
tive  comparison  is  difficult  and  probably  not  worth  while  at  the  present 
stage.  This  model  is  only  one-dimensional  and  therefore  does  not  take  into 


account  the  diffraction  of  the  emitted  radiation  that  will  occur  in  the 


crystal.  This  is  a  rather  serious  shortcoming  of  the  present  model  and  is 
being  rectified. 


IS 


APPENDIX  B 

Das  Gupta  has  recently  begun  an  experiment  where  the  copper  sample  is 
excited  with  16.997  KeV  x-rays  instead  of  8.95  KeV  x-rays.  The  preliminary 
results  of  this  experiment  show  that  the  nonlinear  response  of  the  Cu  K 
intensity  to  the  pump  beam  intensity  is  even  greater.  The  quantity  8.95  KeV 
is  the  energy E .  necessary  to  ionize  the  copper  atom  by  removing  one  electron 
from  the  Is  shell.  The  16.997  KeV  energy  is  the  sum  of  Ei  and  Ej,  =  8.047 
KeV  which  is  the  energy  of  a  Cu  K  photon.  Because  the  pump  beam  is 
incoherent,  the  experiment  with  E  =  16.997  KeV  can  be  described  in  terms  of 
two  decoupled  two-level  systems  with  swept  excitation  as  follows: 


where  E4  -  Ej  =  16.997  KeV,  'fi(u2  -  u)1)=-R(a)4  -  w  3  )  =  8.047  KeV  and 
u 3  -  tdj  =  8.950  KeV.  The  two  level  system  consisting  of  states  2  and  1 
is  the  same  as  was  considered  in  Appendix  A.  States  4  and  3  are  in  the 
continuum.  Although  these  two  states  should  be  represented  by  a  band  of 
states,  we  will  consider  only  a  single  state  in  each  region. 

This  system  of  two  two- level  atoms  is  described  by  the  equations 
of  motion  for  the  density  matrix  elements  and  the  wave  equation 
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P22  =  Mz)  (t  -  ~)  “  (Ya+Y2)P22  -  ft  CV2  IP  1 2  -  Vl2P2l),  (B-l) 


Pll  =  YaP22  “  YlP 1 1  -  ft  (Vl2P21  -  V2lPl2)>  (B.2) 


Yl+Y2+Ya  i 

p2l  =  -  i(t»2  + - 2 - )p21  +  tT  ^P22  "  pl^V21  *  (B  •  3) 


P44  =  A(z)6(t  -  ft)  -  Y4P44  -  ft  (V43P34  -  V34P43),  (B.4) 


P33  =  -  Y3P33  ~  ft  (V34P4 3  -  V43P34)  , 


(B.5) 


P43 


=  -  i(w2  +  Y34)P43  +  fT  ^p44  -  P33)v43 


CB.6) 


P 12  =  P21  j  P34  =  P43  , 


(B.  7) 


( -  -  —  - — )E(z,t)  =  uo  —  P(*.t)  » 
3z2  c2  3t2  at2 


(B.8) 


P(z,t)  =  N  [y2i(p21  +  P12)  +  P43CP43  +  P34)]  , 


(B.9) 
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(B. 12) 


rr 


P33  “  P 3 3 (0) e 


(B. 13) 


where  5  =  t  -  —  and  6(C)  is  the  unit  step  function.  At  time  t  =  0  there 
is  no  population  in  the  state  3  so  p33(0)  =  0  and  P33(t)  =  0.  The 
second  term  in  the  brackets  in  equation  (B.ll)  therefore  becomes 


y2  f  dT'e"Y34(t"t'^(z,t')9(?) 

43  J 


,  -Y45  ft  —5^  T 

=y2  e  9(0  dte  z  «?(z,t-x) 
43 

0 


~Y4?  ft 

-  y?  e  <?(z,t)6(0  dxe 
43  J 


Y4-Y3 


The  remaining  steps  are  identical  with  those  of  Appendix  A  so  they 
will  not  be  repeated.  The  final  solution  is 

(  rz'  f2'  ) 

*2(*',t)  =<?2(0,O  exp/j  r(x)dxF(0  +  rc(x)dxFc(0> f  (B.14) 

where  r(x)  and  F(£)  are  defined  after  Eq.  (A. 13)  and  by  Eq.  (A. 17) 


respectively  and 


r  (x)  =  , 

cv  1  ce0h 


(B. 15) 


Fc(t) 


,r  B.II3.  T 
e'Y4e0(of  dxe  2 


(B. 16) 
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The  only  difference  between  Eq.  (B.14)  and  (A. 17)  (pumping  with  a 
16.997  KeV  instead  of  a  8.95  KeV  beam)  is  the  gain  coefficient  which  is 
increased  by  an  amount 


^Fc(C)j  f(x)dx 


where  A(z)  =  Af(z).  The  k  radiation  still  increases  exponentially  as 

^  1 

the  pumping  rate  A  is  increased. 

This  model  suffers  from  the  same  problems  that  were  listed  in 
Appendix  A.  It  suffers  from  the  further  problem  that  states  4  and  3 
should  be  replaced  by  a  distribution  of  states. 

Comparison  of  the  predictions  of  this  model  and  the  experimental  results 
awaits  completion  of  the  experiment. 
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APPENDIX  C 

In  order  to  comprehend  the  results  of  the  Das  Gupta  experiments,1 
it  is  necessary  to  fully  understand  how  a  crystalline  lattice  affects  the 
processes  that  occur  within  a  single  atom  in  the  crystal.  For  example, 
an  electromagnetic  wave  in  a  crystal  is  more  likely  to  propagate  in  certain 
directions  than  in  others.  This  means  that  any  process  where  the  atom  emits 
radiation  will  be  affected  by  the  crystal.  Of  particular  importance  for  the 
Das  Gupta  experiment*  is  the  effect  of  the  crystal  on  the  spontaneous  and 
stimulated  emission  of  an  atom  in  that  crystal. 

Lee  and  Lee3  have  investigated  the  spontaneous  decay  rate  of  an  excited 

atom  in  a  crystal  made  of  identical  atoms,  resonant  with  the  excited  atom. 

This  is  not  the  case  that  occurs  when  the  excited  atom  is  created  by  the 

removal  of  one  of  the  inner  shell  electrons.  In  order  to  treat  this 

particular  case,  we  have  considered  the  spontaneous  emission  of  an  initially 

excited  two- level  atom  in  a  crystalline  array  of  two-level  atoms.  The 

initially  unexcited  atoms  are  identical  with  each  other  but  not  resonant 

with  the  excited  atom.  The  excited  atom  is  assumed  to  have  a  Bohr  frequency 

too  and  to  located  at  x.  The  remaining  N-l  atoms  have  a  Bohr  frequency 

to  and  are  located  at  positions  x.,  j  =  1,  2,  N-l. 

s  j 

The  equation  of  motion  of  the  probability  ampltiudes  for  the  various 
states  of  the  system  are  given  by 

%  r*  4 

(i  It  "  V  a(t)  =  e  E  A*e  Mt),  (C-l) 

t  k  k 


(i 


3_ 

at 


-  V 


a^(t)  = 


V 


ik-x. 

X(t) 

k 


(C.  2) 
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(i  If  -  wk)bk>(t)  =  r 


.  - i J< * x 

A^c  a(t) 

k 


N- 1 
7  B  e 
j=l  k 


-ik'x. 
3 


(t) 


(C.3) 


with  the  initial  condition  a(0)  =  1,  a.(0)  =  b  (0)  =  0  where  eA  and  eB  are 
•  J  k  k  k 

coupling  constants  proportional  to  dipole  matrix  elements  for  the  initially 

excited  atom  and  unexcited  atoms  respectively,  a(t)  is  the  probability  amplitude 

that  the  initially  excited  atom  is  excited  at  time  t,  all  other  atoms  are  in 

their  ground  state  and  there  are  no  photons  in  the  radiation  field;  b  is  the 

-  t 

probability  amplitude  that  there  is  one  photon  with  eave  vector  k  in  the 

radiation  field  and  all  atoms  are  in  their  ground  state,  and  a^  is  the 

amplitude  for  the  ith  atom  to  be  excited,  no  photons  and  all  other  atoms 

in  their  ground  state.  The  quantity  cA^  is  given  by 

k 


cA  = 

t 


2irhw, 


P/m  +> 


hw, 


(C.4) 


where  V  is  the  volume,  of  the  crystal,  |4>  is  the  ground  state  of  the 
.  .  .  .  -> 

initially  excited  atom  and  P  is  the  momentum  of  the  electron  around  the 
atoms.  A  similar  expression  can  be  written  for  eB^. 

Approximate  solutions  for  Eqs.  (C.l)-(C.S)  can  be  determined  with 
the  aid  of  the  perturbation  theory  of  Lindstcdt  and  Poincard.^  In  this 
method  the  terms  secular  in  t  that  arise  in  conventional  perturbation  theory 
are  used  to  renormalize  the  frequencies  wo,  ii>s,  and  If  the  renormalized 

frequencies  corresponding  to  coQ,  w  ,  and  are  denoted  as  (2, 
respectively  then  new  time  scales  can  be  defined 
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■v  =  Qt , 

(C.5) 

T.  =  ft.t. 

(C.6) 

J  J 

Tk  = 

(C.7) 

so  that  Eqs.  (C.1)-(C.3)  become 

(in  Jj-  -  «  )a(T)  =  e  l  A_*eik'Xby(Tk),  (C.8) 

■f  k  k 

K 


(ift.  -3 - w  )a.(t.)  =  c  T  B  e1  XJb  (t.  )  , 

V  1  dr  sJ  i  £  ?  k 


(C.9) 


(ink  3^  -  “k)bj<Tk>  ■ £ 


A  e 


-ik-x 


N"!  -ith 


a(t)  +  X  BjT1K  Xja  (r .) 
j  =  l  k  3  3 


0 


The  renormalized  frequencies  along  with  the  amplitudes  are  then  expanded 
in  powers  of  e  as 


(1)  ?  (2) 

n  =  U1Q  +  en  +  e2JT  '  +  ...» 

(C. 11) 

(1)  2  (2) 

n.  =  u  +  cQ .  +  c2n.  + 

J  s  j  ; 

(C. 12) 

o  -  +  +  20(2) 

nk  “  wk  G  k  e  ^k  ■*■•••» 

(C. 13) 

(0)  (1)  2  (2) 
a  =  av  '  +  ea^  1  +  e^a'-  J  +  .... 

(C. 14) 

(0)  (1)  2  (2) 
a.  =  a.  +  ea.  +  e2a  '  +  ...» 

J  J  J  j 

(C. 15) 

h>=  bf0>  *  rb™  +  r’h™ 
k  k  k  k 


.10) 


(C.  16 
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When  Eqs.  (C.  1 1) -  (C.  16)  are  substituted  into  P.qs.  (C.8)-(C.]0)  and  the 
coefficients  of  like  powers  of  are  equated  on  the  following  infinite  set 
of  first  order  differential  equations  arises: 


e°  :  (i  j-  -l)a(0)  =  0  , 


l)a<°> 


d 

(lS 


t  k  k 
k 


r  0(n-m)  d  (m)  ,  . 

i  1  fi  7TT  a  C1)  » 


d 

w .  (l  j — 

1  v  dt . 
3 


l)a^  (t  .) 

J  3 


-y  -»- 

*  r  ,•>  lk-X. 

i  Vfn  )(\>  J 
t  k  k 
k 


•  V  r.(n-m)d  (m)  ,  .. 

-i  )  Si)  —  a  (t  .)  , 

n  3  dx .  y 

m=0  J  j 


d  (n),  .  .  (n-1),  .. 

“k(l  57-  -  1)bk  (Tk>  *  V  (T) 

n  k 


N-1  .«►">-  r  1  n 

j‘l  it  3  3 


.  V  i"-"0  d  .fml.  . 

-  ’  ,,  H  ;rr  b  <’d  ■ 


where  n  =  1,2,3, , 


(C. 17) 


(C. 18) 


(C. 19) 


(C. 20) 


(C. 21) 


(C. 22) 
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The  solutions  of  Eqs.  (C.  17) - (C.  18)  are  then 


.(0) 


=  Co 


- 1 1 


-dr 

C.e 

3 


j 


-it 

D  e 


k 


(C. 23) 

CC.24) 

(C. 25) 


where  C,  C.,  D  are  constants  to  be  determined  from  the  boundary  conditions. 

J  t 

Substitition  of  Eqs.  (C.23)  and  (C.24)  into  Eq.  (C.20)  for  the  case  when 
n  =  1  gives: 


w  1 
o 


dt 


(a^ (t)e1T 


_  iOf.)  i V.v 

)  =  ADe  k  e  k  X 

J  L  -V  •> 

k  kk 


-  Cft 


CD 


CC. 26) 


When  this  equation  is  integrated  the  second  term  on  the  right  hand  gives 
rise  to  secular  divergent  behaviour  in  t.  In  order  to  remove  this 
problem  the  solution  is  renormalized  by  choosing  ft 


ft^  =  0. 


Equation  (C.26)  can  now  be  integrated  to  give: 

„(D 


★ 

AD 

i  ft  y  k  £  1Tk  ilCx 

CO  =  l  —  e  e 


mo  J 


ft.  -ft 


Likewise  Eqs.  (C.21)  and  (C.22)  for  the  case  when  n  =  1  can  be  integrated 
to  give 


/ 
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«£»  -  f  >  -  0  , 


* 

B  D  .  t*  -*• 

?  t  IK  ‘  X .  -IT, 
(l)  *  •,  i  r  k  k  ^  j  k 

a.  '(t.)  =  —  )  7T — —  e  Je 

i  v  1  (0  £  ft,  -ft. 

J  J  o  k  k  j 


k  to, 
v  k 


AC  -> 

-ik*x  . 
k  -it 

e  e 


ft-ft 


N-l  C.  -ik-x.  -iT. 
k  J=1  J  k 


In  this  manner,  Eqs.  (C.  20) - (C. 22)  can  be  integrated  for  n  £_2. 

For  each  value  of  n  the  known  quantities  are  suostituted  into  the 
equations  which  are  then  integrated.  Those  terms  that  are  secular  are 
eliminated  from  the  equation  by  setting  their  sum  equal  to  zero.  This 
process  renormalizes  the  frequency. 

Rather  than  write  down  all  the  intervening  steps,  it  is  appropriate 
to  skip  to  the  final  result  because  the  intermediate  steps  are  straight 
forward  but  tedious.  Because  the  object  of  this  exercise  is  to  determine 
the  lifetime  of  an  excited  atom  in  a  nonresonant  crystal,  only  the 
frequency  of  the  solution  of  a(T)  is  important.  This  frequency  to  fourth 
order  in  e  is 


_  N-i  n,n. 

n-n  *  -  —  l  I 

s  u_  £  ..  2 


I A  I  2  I  B 


s  k  j=i  (jysycn-^v 


+ 


N-l  ft 

n  i 


j  %  v 


O'  0),  UK  , 

t  i-i  s  k  k 


*X  C  K 

k  k'  k'  k 


(nk-a)(nk,-n)((ij-n) 


i(it-ic')'  Cx-x.) 
e  } 


4 


(C. 27) 
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where  52^  is  the  renormalized  frequency  that  would  occur  if  the  atom 
were  in  free  space.  In  order  to  simplify  Bq.  (C.27)  the  renormalized 
frequencies  on  the  right  hand  side  can  be  set  equal  to  their  values  when 
the  coupling  constants  are  zero 

fl  -  wQ,  (C.  28) 

Oj  “  «os,  (C.  29) 

°k  *  "k  • 

With  these  approximations  Eq.  (C.27)  becomes 

e  4  l  A_^  | 2  |  | 2 

k  k 

(us-«kHv“k)2 


q  «  -  (N-i)  y  — 

s  i  \ 

k  K 


+ 


1 1 
k  k 


*  * 

A  B  B  B 

t  t'  It  it1 _ 

((VV(V"a)i)(VlV) 


N-l 

l 


3  =  1 


i(it-jt' )  (x-x.) 
e  J 


(C. 31) 


The  summation  over  j  that  occurs  in  the  last  term  of  Eq.  (C.31)  can  now  be 
performed 

N-l  i  (it- it' )  •  (x-x . ) 

l  e  J  “  (N-l)  l  6^  ^  ^  (C.  32) 

j=l  q  k-k',q, 

where  q  is  a  reciprocal  lattice  vector.  The  argument  of  the  delta  function 
in  Eq.  (C.32)  simply  states  the  conditions  for  Bragg  scattering  of  the 
radiation  in  a  crystal. 
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Equations  (C.31)  and  (C.32)  represent  the  current  state  of  these 
calculations.  Future  work  will  include  division  of  Eq.  (C.31)  into  its 
real  and  imaginary  parts  (separate  the  level  shift  from  the  lincwidth) 
and  approximate  evaluation  of  the  sums. 
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The  work  presented  in  this  report  was  performed  by  M.  0.  Scully, 
Professor  of  Physics  and  Optical  Sciences,  University  of  Arizona;  G.  T. 
Moore  and  J.  K.  Mclver,  Research  Physicists,  University  of  Arizona. 

The  part  of  this  work  concerned  with  the  lifetime  of  an  excited  atom 
in  a  crystal  was  presented  in  a  talk  entitled  Spontaneous  Emission  in  a 
Crystalline  Array  at  the  Eleventh  International  Quantum  Electronics 


Conference. 
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